Suzaku observes weak flares from IGR J17391-3021 
representing a common low-activity state in this SFXT 


A. Bodaghee and J. A. Tomsick 

Space Sciences Laboratory, 7 Gauss Way, University of California, Berkeley, CA 94720, USA 

bodaghee0ssl.berkeley.edu 

and 

J. Rodriguez and S. Chaty 

Laboratoire AIM, CEA/IRFU - Universite Paris Diderot - CNRS/INSU, 

CEA DSM/IRFU/SAp, Centre de Saclay, F -9 1191 Gif-sur-Yvette, France 

and 

K. Pottschmidt 

CRESST and NASA Goddard Space Flight Center, Astrophysics Science Division, Code 661, Greenbelt, MD 20771, 

USA 

Center for Space Science and Technology, University of Maryland Baltimore County, 1000 Hilltop Circle, Baltimore, 

MD 21250, USA 
and 

R. Walter 

INTEGRAL Science Data Centre, Universite de Geneve, Chemin dEcogia 16, CH-1290 Versoix, Switzerland 
Observatoire de Geneve, Universite de Geneve, Chemin des Maillettes 51, CH-1290 Sauverny, Switzerland 

and 

P. Romano 

INAF, Istituto di Astrofisica Spaziale e Fisica Cosmica, 

Via U. La Malfa 153, 1-90146 Palermo, Italy 


ABSTRACT 

We present an analysis of a 37-ks observation of the supergiant fast X-ray transient (SFXT) 
IGRJ17391-3021 (=XTEJ1 739-302) gathered with Suzaku. The source evolved from quiescence to 
a low-activity level culminating in three weak flares lasting ~3 ks each in which the peak luminosity 
is only a factor of 5 times that of the pre-flare luminosity. The minimum observed luminosity was 
1.3 x 10 33 ergs -1 (d/2.7 kpc) 2 in the 0.5-1 OkeV range. The weak flares are accompanied by significant 
changes in the spectral parameters including a column density (Ah = (4.1+^) x 10 22 cm -2 ) that is ~2-9 
times the absorption measured during quiescence. Accretion of obscuring clumps of stellar wind mate- 
rial can explain both the small flares and the increase in Ah. Placing this observation in the context of 
the recent Swift monitoring campaign, we find that weak-flaring episodes, or at least epochs of enhanced 
activity just above the quiescent level but well below the moderately bright or high-luminosity outbursts, 
represent more than 60±5% of all observations in the 0.5-10keV energy range making this the most 
common state in the emission behavior of IGR J17391-3021. 

Subject headings: accretion, accretion disks ; gamma-rayl: general ; stars: neutron ; supergiants ; X-rays: binaries ; 
X-rays: individual (IGR J17391-3021 = XTEJ1 739-302) 



1. Introduction 

Supergiant Fast X-ray Transients (SFXTs: Smith 
et al. 2004; Sguera et al. 2005, 2006; Negueruela 
et al. 2006b) form a new subclass within the group 
of high-mass X-ray binaries (HMXBs). Previously, 
HMXBs were catalogued (e.g. Liu et al. 2000, and 
references therein) into one of either two groups: the 
transient sources with B-emission line stellar compan- 
ions (BeXBs); or the persistent, but variable, HMXBs 
with supergiant O or B-type companions (SGXBs). 
The SFXTs share properties from both groups: like 
BeXBs, they are characterized by outbursts of short 
(hours) to long (days) duration whereby the peak in- 
tensity is several orders of magnitude greater than dur- 
ing quiescence; and like SGXBs, they are systems 
in which a compact object (usually a neutron star) 
is paired with a supergiant OB star. Their study is 
therefore important because SFXT s could represent an 
evolutionary link between BeXBs and SGXBs (e.g. 
Negueruela et al. 2008; Chaty 2010; Bodaghee et al. 
2010; Liu etal. 2010). 

XTEJ 1739-302 (Smith et al. 1998) was the first 
hard X-ray transient to display behavior that was iden- 
tified as characteristic of the class of SFXTs (Negueru- 
ela et al. 2006b). It was listed as AX J1 739. 1-3020 
in the ASCA catalog of X-ray sources in the Galac- 
tic Center region (Sakano et al. 2002). Five years 
after its discovery, INTEGRAL detected a transient 
source named IGRJ1 739 1-3021 which, because of 
the positional coincidence and similar emission be- 
havior, was shown to be the hard X-ray (18-50keV) 
counterpart to the RXTE source (Sunyaev et al. 2003; 
Smith et al. 2003). A refined X-ray position was 
obtained with Chandra (Smith et al. 2006), which 
led to the identification of USNO-B 1.0 0596-0585865 
(= 2MASS J17391 155-3020380), a supergiant star of 
spectral type 08 lab located -2.7 kpc away, as the op- 
tical/IR counterpart (Negueruela et al. 2006a; Chaty 
et al. 2008; Rahoui et al. 2008). 

This source has been detected numerous times in 
outburst and in quiescence thanks to monitoring cam- 
paigns with INTEGRAL (Lutovinov et al. 2005; Sguera 
et al. 2005; Tiirler et al. 2007; Blay et al. 2008; Chen- 
evez et al. 2008; Romano et al. 2009b; Ducci et al. 
2010) and with Swift (Romano et al. 2008a, b; Sidoli 
et al. 2008, 2009a, b; Romano et al. 2009a, 2010). Most 
of the outbursts are short (~0.5 h) and weak (~10cps 
in 20-40 keV as seen with INTEGRAL ), while a few 
flares last longer than an hour (sometimes days as seen 


in the 0.5-1 OkeV band with Swift) and are rather in- 
tense (-60 cps or -400 mCrab with INTEGRAL) (Blay 
et al. 2008). The photoelectric absorption is variable 
from outburst to outburst (Smith et al. 2006), reach- 
ing up to (13+ 4 ) X 10 22 cm -2 during the rising portion 
of certain flaring episodes (Sidoli et al. 2009a, b). The 
dynamic range between quiescence and peak outburst 
is - 10 4 (in’t Zand 2005). Based on multiple observa- 
tions of IGRJ1 739 1-3021 gathered over the course of 
a year with Swift , Romano et al. (2009a) showed that 
bright outbursts account for only 3-5% of the entire 
observation set, while the source spends 39±5% of the 
time in an inactive state near quiescence. Quiescence, 
i.e., a flux in the 0.5-1 OkeV range below the Swift de- 
tection limit of 1.6 x 10 -12 erg cm -2 s _1 for a -1 ks ex- 
posure, is a very rare state (Romano et al. 2009a). 

What triggers the outbursts is still an open ques- 
tion. Several mechanisms have been proposed to ex- 
plain the peculiar X-ray emission behavior of SFXTs. 
The outbursts could be due to one or a combination 
of the following effects: the intermittent accretion 
of clumpy wind material by a compact object in a 
larger orbital radius than in classical SGXBs; passage 
of the compact object in an eccentric orbit through 
anisotropic winds emanating from some supergiant 
stars; and the dynamics of wind accretion onto highly- 
magnetized neutron stars including gating effects, 
photo-ionization, transient disks, and Rayleigh-Taylor 
instabilities (e.g. in’t Zand 2005; Walter & Zurita 
Heras 2007; Sidoli et al. 2007; Negueruela et al. 2008; 
Bozzo et al. 2008; Ducci et al. 2009, and references 
therein). 

While the identity of the stellar companion is 
known, the nature of the compact object has yet to 
be conclusively determined. The spectral properties of 
IGR J17391-3021, i.e. its hard power-law shape and 
possible spectral cutoff around 13keV (Sidoli et al. 
2009a), suggest a neutron star primary. However, co- 
herent pulsations or cyclotron absorption lines that 
would firmly establish the presence of a neutron star 
have not been detected thus far (Smith et al. 1998; 
Sguera et al. 2005; Blay et al. 2008). Recently, the or- 
bital period of the system was shown to be 51. 47(2) d 
with hard X-ray outbursts (>20keV) and quiescent 
epochs coexisting at periastron suggesting clumpy, 
anisotropic winds (Drave et al. 2010). 

In February, 2008, Suzaku observed IGR J17391— 3021 
for 37 ks of effective exposure time. The data from this 
observation are described in Section 2. Results from 
timing and spectral analyses are reported in Sections 
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Fig. 1. — Image of IGRJ17391-3021 as captured by the XIS1 detector on Suzaku. Pixels have been grouped in blocks 
of 4 (bin: group4 in ds 9). The circular source and square background selection regions are also shown. Coordinates 
are given as equatorial R.A. and Dec. (J2000) where North is up and East is left. 


3 and 4, respectively. We discuss our findings and 
present our conclusions in Sections 5-6. 

2. Observations 

The Suzaku space telescope features two co-aligned 
instruments working in tandem to observe the hard X- 
ray sky (0.2-600 keV): the X-ray Imaging Spectrome- 
ter (XIS: Koyama et al. 2007) and the Hard X-ray De- 
tector (HXD: Takahashi et al. 2007). Currently, XIS 
consists of two front-illuminated CCD cameras (XIS0 
and XIS3), and one that is back-illuminated (XIS1). 

Observation ID 402066010 (PI: J.A. Tomsick) tar- 
geted IGRJ17391-3021 on February 22-23, 2008 
(MJD 54518.495-54519.396), for 72.271 ks. After 
removing bad events during data reduction (see be- 
low), and observation gaps due to the occultation of 
the satellite by the Earth, the net exposure time that 
remained was 36.546ks. The XlS-nominal pointing 
mode was used with the quarter-window option in or- 
der to gain a higher time resolution. 

Data reduction employed v.6.9 of the HEAsoft soft- 
ware package following the procedures described in 
the Suzaku ABC Guide v.2. The task xi spi converted 
the unfiltered XIS event files to pulse-invariant chan- 
nels, which were then input into xisrepro to re- 
turn cleaned event files. We extracted the source spec- 
trum from each XIS detector using a circular region of 
120"-radius centered on the position given by Chan- 
dra (see Fig. 1). The background was extracted from 
a pair of 240"-wide squares positioned on a source- 
free region of the detector. Based on these source and 


background extraction regions, we generated response 
matrices with xisrmfgen and xissimarfgen. 

To improve photon statistics, we merged the source 
spectral files for the front-illuminated CCDs (XIS0 and 
XIS3), and we did the same for the background. Each 
bin in the XIS spectra contained a minimum of 150 
counts. Source and background light curves were gen- 
erated in 0.5-1 OkeV (XIS) with a time resolution of 
2 s. Timing data from the 3 XIS detectors were com- 
bined into single light curves for the source and for the 
background. 

A tuned (v.2.0) background file and the common 
GTI were used during spectral extraction of HXD- 
PIN. Dead-time corrections were performed with 
hxddt cor. The cosmic X-ray background was mod- 
eled using the flat-field response file from January 29, 
2008, and its contribution was accounted for in the PIN 
background spectrum. Since IGRJ1 739 1-3021 was 
very weak during this observation, PIN recorded few 
high-energy photons. Furthermore, IGRJ1 739 1-3021 
is located close to the Galactic Center so diffuse emis- 
sion from the Galactic Ridge can influence the high- 
energy spectrum of IGRJ17391-3021 (Yuasa et al. 
2008). Applying the best-fitting parameters for the 
ridge emission from Valinia & Marshall (1998) to the 
PIN spectrum, we find that the ridge contributes over 
50% of the PIN photon counts. Hence, the statistics 
being insufficient for timing and spectral analyses, ob- 
servations with PIN are not discussed any further. 
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Table 1 : The different states in the emission behavior of IGRJ1 739 1-3021 as used in the text. 


state 

absorbed flux (0.5-1 0 keV) 
(erg cm -2 s -1 ) 

example time 
(MJD) 

reference 

quiescence 

low 

medium 

high 

< (1 - 2) x 10 -12 
~ (2 - 1 0) X 1 0 -12 
~ (2 - 100) X 10 -11 
> 10- 9 

54518.495-54518.920 

54518.920-54519.396 

54635 

54692 

this work, see also Bozzo et al. (2010) 

this work, see also Sidoli et al. (2008); Romano et al. (2009a, 2010); Bozzo et al. (2010) 
Romano et al. (2009a, 2010) 

Romano et al. (2009a, 2010) 
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Fig. 2.— Light curve of IGR J1 739 1-3021 fromXIS 
(0.5-10keV). The top panel displays the background- 
subtracted light curve from the source, and the bot- 
tom panel shows the background count rate during 
the observation. Also presented are the source light 
curves from soft (S) 0.5-2 keV and hard (H) 2-10keV 
bands, as well as the hardness ratio which is defined 
as (H-S)/(H+S). Each data point collects 240 s worth 
of exposure. After 33 ks of quiescent emission, a 
low state begins at approximately MJD 54518.92 (in- 
dicated by the vertical dashed line). 

3. Timing Analysis 

Figure 2 presents the light curve of IGR J17391 -3021 
in the 0.5-1 0 keV band from XIS. In this energy range, 
the average source count rate (combining XIS 0, 1, 
and 3) for the entire observation is 0.251 ±0.005 cps. 
The best-fitting model to the (absorbed) spectrum of 
the time-averaged Suzaku observation (see Sect. 4) 
yields a coefficient of 2.0 x 10" 11 ergcm" 2 ^" 1 which 
we can use to convert instrumental count rates to 
units of 10" 11 erg cm -2 s" 1 . Thus, the average flux 


is (5.0 ± 0.1) x 10" 12 erg cm" 2 s" 1 . 

Based on Fig. 2, we can identify two distinct emis- 
sion states from IGR J1 739 1-3021 in this observation. 
During the initial 33 ks, IGR J1 739 1-3021 is in a qui- 
escent state in which the average source count rate 
is 0.100±0.009cps ((2.0 ± 0.2) x lO-^ergcnrV 1 ). 
As with previous observations of IGRJ17391-3021 
when it was at or near quiescence (Sidoli et al. 2008; 
Romano et al. 2009a, 2010; Bozzo et al. 2010), our 
Suzaku observation shows that even in this state, 
IGR J17391 -3021 varies on time scales of a few hun- 
dred seconds. 

Around 33 ks after the start of the observation (cor- 
responding to MJD 54518.92), the source enters a low 
state in which the mean count rate of 0.358±0.008cps 
((7.2 ± 0.2) x 10" 12 erg cm" 2 s" 1 ) is a factor 3-4 that of 
the flux in the quiescent state that preceded it. The ob- 
servations prior to and after MJD 545 1 8.92 are hence- 
forth referred to in this analysis as “quiescent” and 
“low” states, respectively (see Table 1). In addition, 
three flares lasting ~2 ks each are revealed in the light 
curve at t ~ 40, 54, and 70 ks. The peak flux of 
1.3±0.2cps ((2.6 ± 0.4) x 10" 11 erg cm" 2 s" 1 ) is an or- 
der of magnitude larger than the flux during the quies- 
cent state, but the maxima are still much lower than the 
typical high-state outbursts expected and seen in such 
objects. We note that the peaks are nearly equidistant 
in time being separated by intervals of approximately 
15ks. 

The light curve was divided into two pieces corre- 
sponding to a soft (S) band between 0.5 and 2keV, 
and a hard band (H) with energies between 2 and 
lOkeV. The hardness ratio can then be defined as 
(H-S)/(H+S) where a negative value indicates a soft 
spectrum, and where a positive value denotes a hard 
spectrum. Figure 2 displays no obvious correlation be- 
tween the hardness ratio and the source intensity. Bin- 
ning the light curve into longer time segments, and 
plotting the hardness ratio versus the source intensity 
shows only a marginal correlation (but within the er- 
rors) between the two parameters at large intensities. 

As part of a monitoring campaign focusing on 
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Fig. 3. — Light curve of IGRJ17391-3021 from Suzaku-X IS, SW/f-XRT, and XMM-Newton. Fluxes are in units of 
10 -11 ergcnT 2 s -1 as observed (i.e. absorbed) in the 0.5-10keV energy band, and time is given in MJD (-54000). 
The top panel displays the background-subtracted light curves from the source during Obs. ID 00030987012 from 
XRT (blue circles, 100-s binning) and from our entire XIS observation (red crosses, 240-s binning). Note the half- 
day gap between the observations. In the bottom panel, the source light curve based on -1 ks observations from the 
XRT monitoring campaign of Romano et al. (2009a) is shown (blue circles) along with their 3cr upper limits for non- 
detections (green downward arrows). The blue upwards arrows designate the maximum flux of high-state outbursts 
detected by XRT that are situated beyond the scale of the graph (from left to right, their values are 125, 12, and 551). 
The average flux from our 37-ks XIS observation is plotted as a single red cross at the intersection of the dotted lines. 
The purple triangle at MJD 54740 corresponds to the average flux in a 31-ks observation with XMM-Newton from 
Bozzo et al. (2010). 
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SFXTs, Swift targeted IGRJ17391-3021 at various 
times throughout 2008 (Romano et al. 2009a). In par- 
ticular, there is a ~900-s long XRT pointing (Obs. ID: 
00030987012) several hours prior to our Suzaku obser- 
vation. The 0.5-10keV Swift-XRT detections and 30 - 
upper limits are shown along with the average count 
rate recorded by XIS in Fig. 3. The instrumental con- 
version coefficient for XRT is 1 .3 x 1 0 -10 erg cm -2 ct -1 
(Romano et al. 2009a). 

4. Spectral Analysis 

The time-averaged spectrum of IGRJ17391-3021 
gathered in the 0.5-1 OkeV range by XIS is presented 
in Fig. 4. We fit an absorbed power law, employ- 
ing the abundances of Wilms et al. (2000) and the 
photo-ionization cross sections of Balucinska-Church 
& McCammon (1992), to the XIS spectra which yields 
T =1.4±0.1 and Ah =(3.6±0.4)xl0 22 cm -2 (errors are 
provided at the 90% confidence level). These are com- 
patible with the values derived by Sidoli et al. (2008) 
for out-of-outburst epochs in a similar energy range 
with Swift-X RT, with an average luminosity that is 
slightly higher at 4.8 X 10 33 ergs -1 compared with 
3.9 X 10 33 ergs -1 for Sidoli et al. (2008). Both stud- 
ies assume a source distance of 2.7 kpc (Rahoui et al. 
2008). The model fits the data well (xl = 0.82 for 66 
degrees of freedom: dot). 

Residuals remain at low energies (< 1 keV) which 
could hint at the presence of soft excess emission. 
We know that this source exhibits a soft excess since 
this emission was detected with XMM-Newton (Bozzo 
et al. 2010). To account for this feature, we added a 
blackbody whose temperature was fixed to the best- 
fitting value from XMM-Newton (kT = 1.2 keV). This 
reduces the xl to 0.55 for 65 dof (F-test probability 
of 5xl0 -5 ), and Ah = (1.71+{|^) x 10 22 cm -2 , but the 
photon index is poorly constrained (T = 0.2+® 3 ). 

We split the spectra into two segments correspond- 
ing to the quiescent and low states described in Sect. 3 
(i.e. before and after MJD 54518.92). Table 2, which 
lists the spectral parameters of IGRJ1 739 1-3021 for 
each epoch, shows that there is a clear distinction in 
column density, photon index, and luminosity. The Ah 
is at least twice as high in the low state ((4.1+21) x 
10 22 cm -2 ) as it is in the quiescent phase ((l.O+jj*) x 
10 22 cm -2 ). The source spectrum is softer (T = 1.5 ± 
0.1) during the low state compared to the quiescent 
state (T = 1.0 ± 0.3). The unabsorbed 0.5-1 OkeV lu- 
minosity rises from 1 .3 when the source is quiescent to 


7.4x1 0 33 erg s 1 when it is in the low state. 

Figures illustrates the incompatibility (at >99% 
confidence) in the absorbing column from quiescence 
to the low state. It is clear from this figure that the col- 
umn density is the main parameter that evolves from 
quiescent to low states. Indeed, for a constant T ~ 1.3, 
the Ah ranges from ~1 to ~4 x 10 22 cm -2 between qui- 
escent and low states, respectively. This is different 
from what was seen by Bozzo et al. (2010) who found 
that the column density remained steady (within the 
statistical errors) during the low state while the T var- 
ied by around 50%. 

To test whether we overestimated the column den- 
sity during the low state by not considering a soft ex- 
cess, we included a blackbody component in the spec- 
tral fits to quiescent and low states, and we fixed its 
temperature to the optimal value from Bozzo et al. 
(2010) (1.2 keV). The absorbing columns remain in- 
consistent (at the 90% confidence level) between qui- 
escence (Ah < 1 .35 x 10 22 cm -2 ) and the low state (Ah 
= (1. 99+027) X 10 22 cm" 2 ). When the blackbody is 
fixed to the lower temperature of 0.13keV provided 
by Bozzo et al. (2010), the column densities are com- 
patible within the large error bars between the quies- 
cent (Ah = (2.5+2q) x 10 22 cm -2 ) and low states (Ah 
= (5.8+ j-j) x 10 22 cm -2 ), but in this case, the result- 
ing blackbody radius is £ 1 00 km which is much larger 
than the expected size of the compact accretor. 

For ease of comparison with previous studies, we fit 
the source spectrum with models used for this source 
by other authors (Sidoli et al. 2008; Romano et al. 
2009a; Bozzo et al. 2010). An absorbed blackbody 
provides a better fit than the absorbed power law (^ = 
0.71 for 66 dof). The parameters from this model 
agree with those describing the blackbody-modeled 
spectrum of IGRJ17391-3021 when it was in the low 
state as seen with Swift (Sidoli et al. 2008; Romano 
et al. 2009a). The results of this model fit to the vari- 
ous states (i.e. time-averaged observation, quiescence, 
and low state) are listed in Table 2. Once again, statis- 
tically significant differences emerge for the column 
density (as well as for the radius of the blackbody 
emission region) for spectra collected in quiescence 
(Ah < 0.22 x 10 22 cm -2 ) and during the low state (Ah 
= (1.4+° 3 ) x 10 22 cm -2 ). 

In the low state, a power law with a partial-covering 
absorber and an additional component to account 
for Galactic absorption (fixed to 1.3 x 10 22 cm -2 
from Kalberla et al. (2005)) yields Ah = (5.1 +J 1 ?) X 
10 22 cm -2 with a covering fraction of 91 ± 4% and 
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Fig. 4. — Background-corrected spectrum of 

IGRJ17391-3021 fit with an absorbed power law. The 
data represent photon counts from XIS1 (red), and 
XISO combined with XIS3 (blue). Each bin collects 
a minimum of 1 50 counts. 


T = 1.76+0 25 (^/dof = 0.56/53). In other words, dur- 
ing the low state that we observed, a large percentage 
of the X-ray emission is shielded by an absorber whose 
TVh is higher than the expected interstellar value. We 
could not constrain the value of the partially-absorbing 
column in the quiescent state. 

A cutoff power law model gives a good fit to the 
time-averaged spectrum (^/dof = 0.57/65). The 
power law index is F = -0.2+[}'g with a cutoff at 
E cut = 3+ 2 keV, and the column density is Ah = 
(2.1 ± 0.3) X 10 22 cm" 2 . These are consistent with 
the parameters listed in Bozzo et al. (2010). The in- 
dex and cutoff energy are compatible between the qui- 
escent (T = -0.1+}'?, E cut > 2keV) and low (T = 
-0.3+jq, E cut = 3+} keV) states. In this case, the 
upper iimit to the column density during quiescence 
(< 1.5 X 10 22 cm -2 ) lies within the 90% confidence in- 
terval of Ah in the low state ((2.2+J 2 ) x 10 22 cm -2 ). 

The spectra were also fit with the Mewe-Kaastra- 
Liedahl (MEKAL) model employing a dual absorber: 
one that affects the system’s soft X-ray emission from 
hot diffuse gas in the supergiant star’s wind (AH Sg ); and 
another that impacts the high-energy power-law emis- 
sion from the compact accretor (Ah x ). The gas tem- 
perature was set to 0.15keV (Bozzo et al. 2010). In 
the quiescent state, we obtain Ansg < 2.2 x 10 22 cm -2 
and A Hx < 2.4 x 10 22 cm -2 for T fixed to the optimal 
value of 1.0 from the absorbed power law fit to the 
quiescent spectrum (^ 2 /dof = 0.54/52). For the spec- 
trum from the low state, Ansg = (5.8 +!}' q) x 10 22 cm -2 


Table 2: Spectral parameters from a power-law (PL), 
a cutoff power law (CP), and a blackbody (BB) model 
fit to XIS 0-3 spectra of IGRJ17391-3021. Fluxes 
(10“ 12 erg cm -2 s -1 ) and luminosities (10 33 erg s -1 ) are 
given as unabsorbed in the 0.5-10 keV range for a 
distance of 2.7 kpc. Quiescence and the low state re- 
fer to pre and post-MJD 545 1 8.92 spectra, respectively 
(see Fig. 2). The Ah is given in units of 10 22 cm -2 , 
the cutoff energies and blackbody temperatures are in 
keV, and the blackbody radii are in meters. Errors are 
quoted at 90% confidence. 


model 

parameter 

time-averaged 

quiescence 

low state 

PL 

AT h 

3.6 ±0.4 


4 1+0.5 
’ -0.4 


r 

1.4 ±0.1 

1.0 ±0.3 

1.5 ±0.1 


F 

5.5 

1.5 

8.5 


L 

4.8 

1.3 

7.4 


^/dof 

0.82/66 

0.68/8 

0.73/54 

CP 

AT h 

2.1 ± 0.3 

< 1.5 

2 2 40,7 

_03^ 

U,J T 1.0 

3!! 


r 

A"cut 

3!? 

-0.1 ±1.3 
>2 


F 

4.0 

1.3 

5.7 


L 

3.5 

1.1 

4.9 


^J/dof 

0.57/65 

0.64/7 

0.44/53 

BB 

AT h 

1.2 ±0.2 

< 0.22 

1 440.3 
, 49*57 
88!? 


kTbb 

Rbb 

1 40+O.O6 
J -0.05 

72 ±5 

1 C7+0.16 

37 )j 


F 

3.4 

1.2 

5.1 


L 

3.0 

1.0 

4.4 


*$/dof 

0.71/66 

0.73/8 

0.51/54 


Table 3: Same as Table 2, but for key parameters from 
rate-selected spectra from epochs of low count rate 
(LCR), medium count rate (MCR), and high count rate 
(HCR). 


model 

parameter 

LCR 

MCR 

HCR 


count rate range (cps) 

<0.1 

0. 1-0.3 

>0.3 


average count rate (cps) 

0.05 

0.19 

0.45 


exposure time (ks) 

22 

11 

4 

PL 

AT h 

<0.9 

4 1+1.0 
■ -0.9 

4 3 +1 -° 
^ -0.8 


r 

^-ol 

1.7 ±0.3 

1.4 ±0.2 


*5/dof 

1.2/11 

0.83/34 

0.50/38 

CP 

AT h 

<0.2 

2 l +u 
’ -0,5 

2.5t| 


r 

-0.8 ± 0.2 

-0 i +0 - 4 
-0.2 

-0.3 ± 0.2 


^cut 

3 (fixed) 

3 (fixed) 

3 (fixed) 


rfl dof 

1.2/10 

0.71/34 

0.27/38 

BB 

AT h 

<0.1 

1 2 +0<5 
1 ’ -0.4 

1 2 +0,5 
1 ' -0.4 


kTyb 

'Ki 

1.4 ±0.1 

1.4± 0.1 


X $! dof 

1.6/10 

0.72/34 

0.72/34 


and A Hx < 1.4 x 10 22 cm -2 for T = 1.8 ± 0.2 (^ 2 /dof 
= 0.55/65). This suggests that not only is there mate- 
rial around the binary system absorbing X-rays, there 
is additional matter located near the compact object. 
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Fig. 5. — Parameter space of Nu and T derived 
from absorbed power laws fit to the spectra of 
IGRJ1 739 1-3021 during quiescent and low states. 

The best-fitting parameters are represented by crosses 
and the contours denote confidence levels of 68%, 

90%, and 99%. 

We sought to determine how the spectral parame- 
ters varied with the source intensity. So we created 
spectra from time intervals in which the count rates 
(CR) were <0.1 cps (“low CR” or LCR), between 0.1 
and 0.3 cps (“medium CR” or MCR), and >0.3 cps 
(“high CR” or HCR). Then, each rate-resolved spec- 
trum was fit with an absorbed power law, an absorbed 
cutoff power law with fixed to our best-fitting 
value (3 keV), and an absorbed blackbody model. Ta- 
ble 3 lists the results of the spectral fits. Regardless of 
the spectral model, the 7V H is significantly smaller dur- 
ing the LCR epochs. The MCR and HCR epochs have 
spectra that are consistent, and both are statistically in- 
compatible with the spectrum from LCR epochs. 

5. Discussion 

5.1. Weak Flares and Low Activity 

Our analysis of a Suzaku observation of IGR J1 739 1-3021 
revealed the source evolving from quiescence to a 
low-activity state. Yet even at this low-flux level, 
the emission from IGR J17391 -3021 is character- 
ized by variability on time scales of hundreds to 
thousands of seconds indicating an accretion that 
is continuous and dynamic. The average luminos- 
ity of IGRJ17391-3021 during the entire Suzaku 
observation is 4.8 x 10 33 ergs -1 (unabsorbed) and 
3.1 x 10 33 ergs -1 (absorbed) in the 0.5-10keV range, 
and assuming a source distance of 2.7 kpc (Rahoui 


et al. 2008). 

During the quiescent state, which consists of the 
first ~33 ks of the observation, the luminosity in the 
0.5-10keV band is only 1.1 x 10 33 erg s -1 as observed, 
and 1.3 x 10 33 ergs -1 when corrected for absorption. 
This is equivalent to the lowest luminosity detected 
by ASCA (Sakano et al. 2002). In our flux-limited 
spectrum (LCR), the unabsorbed luminosity is 8.7 x 
10 32 erg s -1 , i.e., twice the minimum values measured 
with XMM-Newton (Bozzo et al. 2010) and Swift (Ro- 
mano et al. 2010). 

Concluding our observation is a ~39-ks epoch of 
enhanced activity including three flares lasting 1-3 ks 
each where the peak luminosities are only a factor 
of ~5 greater than the quiescent baseline emission 
that preceded it. Clearly, these flares are significantly 
weaker than the typical high-luminosity outbursts of 
SFXTs in which the dynamic range is a few orders of 
magnitude. These small spikes are reminiscent of the 
short, low-intensity flares seen in this source by Bozzo 
et al. (2010) with XMM-Newton. 

Separated by regular intervals of approximately 
15ks, the flares appear to be quasi-periodic. This 
timescale is well short of the established orbital pe- 
riod of 51 days (Drave et al. 2010), and it is prob- 
ably too long for a pulsation period if the compact 
object is, indeed, a neutron star. Ducci et al. (2010) 
detected quasi-periodic outbursts of similar frequency 
from this source in INTEGRAL observations, but for 
longer flares (between 5 and lOks) of much higher 
intensity (L x ~ 10 36 ergs -1 ). They ascribe the quasi- 
periodic flaring to the formation and dissipation of 
an unstable accretion disk which is the mechanism fa- 
vored by Smith et al. (2006) to explain bright outbursts 
in this source. The pulsation period of the suspected 
neutron star has never been measured, despite numer- 
ous observations in and out of outburst. If the pulse 
period of the suspected NS in this system were known, 
then a change in the sign of the period derivative could 
indicate a dissipation and reformation of a disk (Taam 
et al. 1988). This 15-ks frequency could possibly cor- 
respond to a periodicity but we would require further 
observations to confirm or reject this possibility. 

5.2. Spectral Evolution 

The evolution of the emission from quiescence to 
the low state is accompanied by changes in the X-ray 
spectral parameters. The TVh is higher and the spec- 
tral slope is steeper during the low state than during 
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quiescence. Similarly, the high count-rate resolved 
spectrum showed larger TVh and photon index than the 
low count-rate resolved spectrum. The main parame- 
ter that evolves is the column density. This is different 
from the weak flares seen by Bozzo et al. (2010) who 
found that as the average count rate of the spectrum in- 
creased, T became smaller (i.e. flatter, harder spectra, 
and more positive hardness ratios) while the Nu did 
not change significantly. 

Differences in the spectral properties of the flares 
caught by Suzaku and the ones seen with XMM- 
Newton could be due to unequal geometric configu- 
rations of the accretion stream between the observa- 
tions. The capture by the NS of a clump of wind mate- 
rial will lead to an increase in the mass-accretion rate 
which both instruments detect as a spike in the X-ray 
luminosity. A partial-covering model fit to the X-ray 
spectra from the weak flares observed by Bozzo et al. 
(2010) suggests that up to 40% of the X-ray emission 
is not absorbed by the clump itself, but rather by the 
interstellar medium. The same model fit to our low- 
state spectrum reveals an absorber that shields a much 
larger fraction of the X-ray source (less than 13% of 
the emission is unblocked). 

On the other hand, if prior to being accreted, the 
captured material crosses the line of sight between 
the observer and the NS, the flare will be accom- 
panied by an increase in the column density, which 
is precisely what we observe (see also Rampy et al. 
2009). The expected interstellar value along the line 
of sight is (1.2-1.3)xl0 22 cnT 2 (Dickey & Lockman 
1990; Kalberla et al. 2005), consistent with what we 
measured during quiescence. This represents an up- 
per limit in the direction of the source since it results 
from the integration of absorbing material through the 
entire Galaxy, while IGRJ1 739 1-3021 is located only 
~3 kpc away. The from a power law fit to the spec- 
trum of IGRJ17391-3021 during its low state is 2-A 
times the interstellar value and 2-9 times the quies- 
cent value. This suggests that the source is not strongly 
absorbed (intrinsically) during this observation except 
when the NS accretes a clump passing along our line of 
sight. Material in a clump located near the NS will be 
photo-ionized so the true column density is expected 
to be higher than our measured value which assumes 
cold neutral gas. Our largest column density remains 
an order of magnitude below the maximum TVh mea- 
sured for this source, (32 ± 3) x 10 22 cm -2 (Smith et al. 
2006), although their result relied on a spectrum from 
RXTE that did not extend below 2.5 keV. 


If the origin of the hard X-ray emission is a black- 
body, then this radius expands from 37^ m during 
quiescence to 88 m when the source is in the low 
state, in agreement with Romano et al. (2009a) who 
showed the blackbody radius growing with intensity. It 
is tempting to associate this small blackbody-emitting 
region with the polar caps of a magnetized accret- 
ing neutron star (i^b ^ O.IRns* Hickox et al. 2004). 
However, Bozzo et al. (2010) caution that models de- 
scribing spherically-symmetric accretion onto magne- 
tized neutron stars predict an inverse relation between 
the size of the hot spot and the square root of the X- 
ray luminosity, at least for systems with an average 
L x < 10 35 ergs" ! (White et al. 1983, and references 
therein). 

5.3. The Ubiquity of Low Activity 

Our Suzaku observation of IGRJ17391-3021, 
when placed in the context of the long-term monitoring 
campaign of this source by SvW/f-XRT (Romano et al. 
2009a), shows that the average flux that we recorded 
is consistent with many of the XRT detections at a low 
flux level (F ~ (2 - 10) X 10“ 12 erg cm -2 s -1 ), and is 
on par with a few of the 3 cr upper limits to the count 
rates of non-detections as measured by XRT (Fig. 3). 

Weak flares that are just above quiescence, and 
whose peak fluxes are at the level of the faint de- 
tections in the Swift monitoring data, compels us 
to propose four states in the emission behavior of 
IGRJ1 739 1-3021: the extremely low but variable 
emission of quiescence (average F a b s < (1 - 2) x 
10" 12 erg cm -2 s" 1 in 0.5-1 OkeV); a low-activity state 
that potentially includes weak flares such as those seen 
here and in Bozzo et al. (2010) (F ~ (2 - 10) x 
10“ 12 ergcm“ 2 s -1 ); a medium state with flares of 
larger amplitude (F ~ (2 - 10) x 10" 11 erg cm" 2 s" 1 ); 
and finally, a high state in which the outbursts repre- 
sent the maximum source luminosity (F £ 1 0" 9 erg cm" 2 s" 1 ). 
The boundaries are not strict since the peak fluxes 
in the low state extend up to 3 x 10" 11 erg cm" 2 s" 1 , 
but they provide a useful classification for the various 
emission states in this source. 

If we assume that the Swift detections when the 
count rates were sufficiently low (i.e. F ~ (2 - 10) x 
10" 12 erg cm" 2 s" 1 ) were due to low states that may or 
may not include weak flares, then we can define a low- 
state duty cycle (DC) as follows: 

DC =60 ±5% (1) 
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Here, £7j ow is the sum of exposure times of all ob- 
servations in which IGRJ1 739 1-3021 is detected at 
an average flux within the range of F ~ (2 - 10) x 
10 - 12 ergcm - 2 s -1 , including the Suzaku and XMM- 
Newton observations, while Zr to t denotes the total ex- 
posure from all observations. As expected, the duty 
cycle is not far from the value of 61 ±5% provided in 
Romano et al. (2009a, 2010). 

Low-activity states in which the flux is above the 
“typical” quiescent level, but well below the medium 
and bright flaring episodes, are therefore quite com- 
mon in this source. In the field of IGRJ17391-3021, 
the limiting flux for XRT is 1.6 X 10 -12 and 1.5 X 
10 -13 ergcnT 2 s -1 for exposure times of 1 and 30ks, 
respectively (Romano et al. 2009a). Therefore, most 
Swift upper limits would be detections if the exposures 
had been equivalent to those of the Suzaku and XMM- 
Newton observations suggesting that the duty cycle of 
the low state is much higher. 

According to the source ephemeris (Drave et al. 
2010 ), the midpoint of our observation coincides with 
an orbital phase of 0.37 ^-q 2 whereas the XMM-Newton 
observation took place at phase 0.8. This is another in- 
dication that the low state is common since it is not 
confined to a specific part of the orbit. 

Despite the large duty cycle, weak flares lasting 
a few hundred to a few thousand seconds can be 
missed by short monitoring observations. Because 
their peak flux is so low, detecting a low state re- 
quires prolonged observations where the long integra- 
tion time guarantees that we will achieve the sensi- 
tivity necessary to detect such low levels of activity. 
Given that IGRJ17391-3021 is a prototypical mem- 
ber of its class, it would not be surprising to find 
weak flares from other SFXTs as telescopes continue 
to devote observing time to such objects. In fact, 
recent Suzaku and XMM-Newton observations of the 
SFXT IGRJ08408-4503 included similar low states 
that were punctuated by weak flares (Bozzo et al. 20 1 0; 
Sidoli et al. 2010). 

6. Summary & Conclusions 

The main results from our Suzaku observation of 
the SFXT IGRJ17391-3021 can be summarized as 
follows: 

- For most of the observation, IGRJ17391-3021 
is in what we believe represents the quiescent 
state. The unabsorbed luminosity in the 0.5- 


lOkeV range is 1.3 x 10 33 ergs _1 , and as low 
as 8.7 x 10 32 ergs _1 in the flux-limited spec- 
trum. The spectrum from this epoch is well de- 
scribed by an absorbed, hard power law with Ah 
= (1.0*};*) x 10 22 cm " 2 and T = 1.0 ± 0.3, or by 
an absorbed thermal blackbody with a tempera- 
ture of 1 .57+°;|* keV. 

- Around 33 ks into the observation (MJD 545 1 8.92) 
the source became slightly more active (“low 
state”) culminating in a series of weak flares 
each lasting ~3ks whose peak luminosity is 
only a factor of ~5 greater than the quiescent 
emission that preceded it. During this low state, 
the Ah reached (4.1 ^ 4 ) x 10 22 cm -2 which is a 
factor of ~2-\ times that of the expected inter- 
stellar value, and which is ~2-9 times the Ah 
measured during quiescence. The accretion of 
small clumps of stellar wind material can be re- 
sponsible for the weak flares, while the increase 
in Ah can be explained by the passage of the 
clump in the line-of-sight to the X-ray source 
prior to accretion. 

- The average luminosity of the entire observation 
is4.8xl0 33 ergs _1 . When placed in the full light 
curve from the Swift-XRT monitoring campaign 
by Romano et al. (2009a), we find that this lu- 
minosity is on par with the low-intensity XRT 
detections (as well as a few of the 3 cr upper lim- 
its from non-detections). Furthermore, we sus- 
pect that many of the faintest detections from the 
XRT contain similar low states (potentially in- 
cluding weak flares) that are just above the qui- 
escent flux. We estimate that the duty cycle of 
the low state is 60±5%, which makes it the most 
common state for this source. 

As a prototypical member of the class of SFXTs, 
IGRJ17391-3021 holds valuable clues to the accre- 
tion processes of these intriguing transients. Contin- 
ued monitoring of this source by sensitive instruments 
such as Suzaku should fuel the ongoing debate con- 
cerning the mechanisms that are responsible for their 
peculiar emission behavior. 
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